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ABSTRACT 

We report the results of a sensitive survey of young planetary nebulae in the 
CO J = 2 — 1 line that significantly increases the available data on warm, dense, 
molecular gas in the early phases of planetary nebula formation. The observa- 
tions were made using the IRAM 30 m telescope with the 3x3 pixel Heterodyne 
Receiver Array (HERA). The array provides an effective means of discriminating 
the CO emission of planetary nebulae in the galactic plane from contaminating 
emission of interstellar clouds along the line of sight. 110 planetary nebulae were 
observed in the survey and 40 were detected. The results increase the number 
of young planetary nebulae with known CO emission by approximately a factor 
of two. The CO spectra yield radial velocities for the detected nebulae, about 
half of which have uncertain or no velocity measurements at optical wavelengths. 
The CO profiles range from parabolic to double-peaked, tracing the evolution 
of structure in the molecular gas. The line widths are significantly larger than 
on the Asymptotic Giant Branch, and many of the lines show extended wings, 
which probably result from the effects on the envelopes of high velocity jets. 
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1. Introduction 

Observations of molecular gas have played a key role in developing our current under- 
standing of stellar evolution from the Asymptotic Giant Branch (AGB) through the planetary 
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nebula (PN) phase. Stars evolving on the upper AGB undergo an increasing rate of mass 
loss which forms a dense, slowly expanding circumstellar envelope of dust and molecular gas. 
The mass loss eventually terminates the evolution on the AGB, and when the surface tem- 
perature of the star becomes high enough to emit a significant flux of ultraviolet radiation 
(T > 30,000 K), it begins to photo-ionize the circumstellar material and form an optically 
visible PN. 

Observations of neutral gas around bona fide PNe have been crucial in the development 
of this general scenario, and numerous atomic and molecular species have been observed. 
These include Hi, Cl, Ol, and other neutral atoms (e.g., Rodriguez, Goss, & Williams 2002; 
Bachiller et al. 1994; Liu et al. 2001; Dinerstein, Sneden, & Uglum 1995), as well as H2 (e.g., 
Kastner et al. 1996), OH (e.g., Zijlstra et al. 2001), CO (e.g., Huggins et al. 1996), and other, 
rarer molecular species such as HCO + and CN (e.g., Bachiller et al. 1997). Millimeter CO 
emission has proved to be the most generally useful probe of the warm molecular component 
of the neutral gas for the same reasons that it is a standard probe of the precursor AGB 
envelopes: CO is relatively abundant and the millimeter transitions are easily excited. It can 
therefore be used to study the mass, structure, and kinematics of the bulk of the molecular 
gas. 

Most of the known CO envelopes around PNe were discovered in dedicated CO surveys 
carried out by Huggins & Healy (1989) using the NRAO 12 m telescope, and by Huggins et 
al. (1996) using the SEST 15 m and IRAM 30 m telescopes. The compilation in Table 3 of 
Huggins et al. (1996) lists 44 PNe detected in CO. Several of these have provided examples 
for detailed study at high angular resolution, and interpretation of the ensemble of the CO 
observations forms an integral part of the evidence for the evolutionary scenario outlined 
above. 

Several recent lines of investigation, including optical imaging of PNe (e.g., Sahai & 
Trauger 1998; Lopez 2003), kinematic studies of proto-PNe (e.g., Bujarrabal et al. 2001), and 
observations of individual molecular envelopes (e.g., He 3-1475, Huggins et al. 2004; AFGL 
618, Cox et al. 2003), have focused attention on the complex geometry of PN formation. 
The underlying mechanisms are not well understood, but they likely involve effects that 
influence the survival of, or affect the characteristics of, the molecular gas found in young 
PNe. Accordingly we have undertaken a deeper search for CO, especially in young PNe, to 
significantly extend the sample where the characteristics of the molecular gas are known. 

Based on considerations of line opacity, atmospheric transparency, telescope beam size, 
and receiver sensitivity, the 230 GHz CO J = 2 — 1 line is the optimal transition for CO 
searches of PN envelopes, where the emission is distributed on small angular scales and is 
optically thin, or partially thin, in most cases. In the time since our previous survey work, 
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an important increase in capability at 230 GHz has been the development of multiple beam 
receivers. PNe are distributed close to the galactic plane, and a sensitive search in CO is 
inevitably contaminated by interstellar CO emission. The use of a multiple beam array on 
the IRAM 30 m telescope for the present survey observations proves to be a very efficient 
way to address this issue. 

The plan of the rest of this paper is as follows. Section 2 describes the sample selection 
and the observations; Sections 3 and 4 present the results of the survey, in the form of a table, 
spectra, and comments on individual objects; and Section 5 gives gives a brief discussion. A 
detailed analysis of the results will given in a separate paper. 

2. The Survey 

2.1. Survey Sample 

The survey PNe were taken from the Strasbourg-ESO Catalog of Galactic PNe (Acker 
et al. 1992) and were chosen to include a major component of relatively young PNe with a 
broad range of spectroscopic characteristics. The PNe were selected to have small angular 
diameters (< the telescope beam size whose FWHM is 12"), relatively high fluxes in Ha 
(with a few exceptions > 10~ 13 erg cm" 2 s _1 , uncorrected for reddening), and consequently 
moderate to high surface brightness in the Balmer lines and the radio continuum. The 
PNe lie close to the galactic plane, and the practical declination limit set by the telescope 
location is ~ —25°. In addition, the PNe were distributed among different spectroscopic 
groups, based on line ratios given in the PN catalog, or updated from the literature: 1) 
[Om] A5007/H/3 < 6.0 ; 2) [Om] A5007/H/3 > 6.0 and [Nil] A 6584/Ha > 1.0; 3) [Om] 
A5007/H/3 > 6.0 and [Nil] A 6584/Ha < 1.0. These line ratios are not very sensitive to 
the reddening, but are subject to uncertainties related to the slit position and size relative 
to the nebula. A fourth group of miscellaneous PNe (see below) was also observed as time 
permitted in the program. 

The angular sizes of relatively young PNe are expected to be quite small at distances 
outside our local neighborhood ~ 1 kpc, which has been explored in CO by earlier obser- 
vations. For example, the well known PN NGC 6720 (the Ring Nebula) which is used as 
an intermediate benchmark below, has an angular diameter of 76" at a distance of 450 pc 
(O'Dell et al. 2002), and would have an angular diameter of < 11" at distances > 3 kpc. 
More compact PNe would have even smaller angular diameters, and such objects are well 
suited to line searches made with single observations centered on the nebula. The average 
surface brightness of the PNe in the Balmer lines or the radio continuum is independent 
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of distance, and is a useful guide to the degree of evolution of the nebulae. For discussion 
here we use the surface brightness Tb at a frequency of 1.4 GHz, using the fluxes measured 
directly at 1.4 GHz from Condon & Kaplan (1998) or the equivalent from the (de-reddened) 
Balmer lines, assuming optically thin emission. The emission at 1.4 GHz becomes optically 
thick before that at higher radio frequencies, but high opacity affects relatively few of the 
PNe and this does not substantially affect our overall discussion. The distribution of Tb for 
the whole PN sample ranges from 4 K to > 10, 000 K, with 1st, 2nd and 3rd quartiles at at 
110 K, 800 K, and 4,200 K, respectively. 

It is well known that the surface brightness of PNe roughly correlates with nebula size 
(e.g., Phillips 2002, and references therein), and is therefore related to the dynamical age of 
the nebulae. We do not discuss this in great detail but note as benchmarks that T B = 3, 500 K 
for the compact PN BD+30°3639 and T B = 65 K for the more evolved PN NGC 6720, and 
their expansion ages based on robust proper motions measured with HST are 800 yr and 
1,500 yr, respectively (Li et al. 2002; O'Dell et al. 2002). 32% of the survey PNe have 
T B > 3, 500 K and 84% have T B > 65 K. Thus the bulk of the observed PNe, with the 
exception of those with the lowest values of Tb, are expected to be younger (and in many 
cases much younger) than ~ 2000 yr. Our selection criteria strongly discriminate against 
large, evolved, low surface brightness PNe. Some of these are known to retain significant CO 
envelopes, as exemplified by the Helix Nebula (Young et al. 1999), Tb ~ 1 K, age ~ 10 4 yr, 
but the objective of the present work is to study the earlier evolution of the envelopes. 

It is known from our previous observations (Huggins et al. 1996) that CO envelopes 
of PNe show large differences which are likely related to the mass of the central stars, 
and possibly additional factors as well. The distribution of the survey PNe over different 
spectroscopic groups is intended to provide information to investigate this further. Of course 
the relation between age and nebula excitation is a complex one: young PNe can have a range 
of stellar temperatures because the speed of the temperature increase of the star depends 
strongly on the stellar mass (e.g., Gesicki & Zijlstra 2000). The nebula excitation is therefore 
not a direct measure of the nebula age, but alone, or in combination with the brightness 
temperature may cast light on, or serve as a predictor of the CO emission. 

Based on the [O ill] A5007/H/3 line ratio the first group (40 objects) selects low excitation 
PNe (excitation classes < 2 on the scale of Webster 1975) with central star temperatures 
< 45,000 K (e.g. Kaler & Jacoby 1991). The second group (31 objects) selects higher 
excitation PNe with strongly enhanced [N II] emission and is likely to include the higher 
mass progenitor stars which pass rapidly through the low excitation phase; the third group 
(28 objects) acts as a control. A small, additional fourth group (11 objects) of miscellaneous 
PNe with incomplete spectroscopy, or with larger angular diameters and typically in the 
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lowest quartile of Tb, was also observed as time in the program allowed. 

2.2. Observations 

The survey observations were made in the CO J = 2 — 1 line at 230.538 GHz using the 
IRAM 30 m telescope at Pico Veleta, Spain, in March 2002 and May 2003. 

The receiver system used was the multibeam Heterodyne Receiver Array (HERA), which 
consists of nine receivers arranged in a regular 3x3 grid with a spacing on the sky of 24". 
The beam size of each "pixel" at 230 GHz is 12" (FWHM). For a gaussian source of size 12" 
(FWHM) centered in the array, a negligible fraction (< 1%) of the signal in the central pixel 
is recorded in the outer pixels. The backend used was the VESPA autocorrelator, which was 
configured to provide 225 x 1.63 kms" 1 channels for each of the nine receivers during the 
2002 observing run, and 409 x 1.63 kms -1 channels during the 2003 observing run. A wide 
spectral coverage is essential because the CO lines of PN envelopes often have broad wings 
and the systemic radial velocities of many of the PNe have not previously been measured at 
optical wavelengths (see Section 5). 

The telescope pointing was regularly monitored during the observations using reference 
sources, and is generally expected to be accurate within ~ 2", although it is somewhat 
more uncertain during periods of high wind or rapidly changing weather conditions. The 
observations were made by switching in azimuth by 4' using the wobbler. The spectra were 
calibrated using the chopper wheel technique, and the intensities are reported here as main 
beam temperatures, using a beam efficiency of 0.51. The velocity scale used for the spectra 
is heliocentric. 

The coordinates of the PNe were obtained from a compilation of optical positions (Acker 
ct al. 1992) and accurate radio positions measured with the VLA (Zijlstra, Pottasch, & 
Bignell 1989; Aaquist & Kwok 1990; Condon & Kaplan 1998). In almost all cases the 
coordinates are known to within ~ 1", and in the least favorable cases to within a few 
arc seconds. The coordinates are therefore not a significant source of uncertainty for the 
observations. Nearly all the PNe lie at low galactic latitude (6), with a mean value of |6| of 
3?8. 

Each PN was observed centered in the array. In all cases where the envelope was detected 
there was no significant envelope emission in the surrounding pixels. Interstellar (IS) lines 
were seen along most lines of sight, and they appear as emission (if they occur in the on 
position) or as dips (if they occur in the off position). In some cases the lines are narrow, 
isolated features and there is little confusion with the PN emission. In other cases, where 
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there are multiple components which may be merged, comparison of all the spectra from 
the array is very efficient in discriminating the IS from the PN components. The spectra 
presented in the next section have been obtained by subtracting the average of the outer 
pixels from the central pixel. This has little effect on the S/N, and is helpful in reducing the 
strength of the IS components, especially in cases where they are fairly uniform or have a 
smooth gradient across the field of the array. 

3. Results 

110 PNe were observed in the survey and 40 were detected. The results of the obser- 
vations are summarized in Table 1. Column (1) gives the PN name in galactic notation 
(ordered by galactic longitude); column (2) the common name; column (3) the rms noise 
level of the observations; columns (4)-(7) the parameters of the CO emission in the cases 
where the PN is detected (see below); column (8) the heliocentric velocity interval over which 
there is significant interstellar contamination; and column (9) comments. 

The spectra of the PNe in which CO emission is detected or tentatively detected are 
presented in Fig. 1. It can be seen that the spectra exhibit a large range in S/N ratio 
and a variety of line profiles. The parameters of the emission given in Table 1 refer to the 
integrated line intensity (J), the heliocentric systemic velocity in CO (V ), the line width 
at zero intensity (AV), and the peak line temperature (T p ). In the spectra where the S/N 
ratio is sufficiently high that the profile shape can be discerned, / is obtained by direct 
integration of the profile, V Q is estimated where possible from the symmetry of the line (e.g., 
the average velocity of the peaks of a double-peaked profile), and AV and T p are determined 
by inspection; because the intrinsic profiles are unknown, larger values of AV from broad 
wings at or below the noise level of the spectra cannot be ruled out. In the spectra with low 
S/N where the profile shape cannot be discerned, the parameters are obtained by a Gaussian 
fit to the line; the line width given in Table 1 is the FWHM, and it is given as a lower limit 
to AV. 



4. Comments on Individual Spectra 

Here we comment on the individual spectra in Fig. 1, point out IS features that appear 
in the profiles, and compare with available optical radial velocities (K P t)- These are taken 
from the compilation by Durand, Acker, & Zijlstra (1998), except where explicitly cited. 
Previous CO observations at other telescopes are also noted. 



-7- 



002.6+04.2 (Th 3-27) Tentative 6-a detection. V Q consistent with V opt (-135.3 ± 
10.1 kms -1 ). The narrow CO spike at V = —13 kms -1 is IS. 

003.8+05.3 (H 2-15) Tentative 5-a detection. V Q consistent with V opt (-63.2+3.7 kms" 1 ). 

008.3—01.1 (M 1-40) Roughly flat-topped profile with severe IS contamination; the 
emission at V > kms -1 is probably all IS. V Q consistent with V opt (—34.4 ± 2.0 kms -1 ). 

008.6— 07.0 (He 2-406) Double-peaked profile. V Q consistent with V opt (+28.2+5.2 kms" 1 ). 

010.1+00.7 (NGC 6537) Strong, triangular profile, with IS contamination on the red 
wing. V Q of peak differs from V opt (—16.9 ± 1.9 kms -1 ) by 13 kms -1 , suggesting some 
structural asymmetry. In deep images the nebula is much more extended than the 10" 
listed in the Strasbourg-ESO catalog, but CO is seen only in the central pixel, and likely 
arises in the dense inner torus whose diameter in HST images is ~ 5" (Huggins & Manley 
2005). A broad CO(1-0) feature at —8 kms -1 (+5 kms -1 LSR) reported by Zhang et al. 
(2000) is significantly different from the spectrum reported here and is likely at least partly 
interstellar. 

011.9+04.2 (M 1-32) Tentative 3-a detection. V Q consistent with V opt (-90.6+7.6 kms -1 ). 

016.0- 04.3 (M 1-54) V of the CO peak is offset by ~ 15 kms -1 from V opt (-47 .4 ± 
4.2 kms -1 ), so the tentative, second CO feature at V = —87 kms -1 may be part of a broad, 
double-peaked profile. 

019.9+00.9 (M 3-53) Probably a flat-topped profile, but strong IS (spike/dip) contami- 
nation of the blue wing. V Q consistent with uncertain V opt (+20.6 ± 14.0 kms -1 ). 

021.1- 05.9 (M 1-63) Double-peaked profile. V Q consistent with V opt (+8.5+10.3 kms -1 ). 

021.7— 00.6 (M 3-55) Strong, parabolic profile, with severe IS (spike/dip) contamination; 
emission at V > 40 kms -1 is probably all IS. The broad wing to the blue appears to belong 
to the PN. V Q consistent with uncertain V opt (+9.7 ± 15.0 kms -1 ). 

021.8— 00.4 (M 3-28) Strong, double peaked profile with IS (spike) contamination on red 
wing and at higher velocities. A CO spectrum obtained with the SEST 15 m telescope was 
previously reported by Gomez, Rodriguez, & Garay (1992). V Q consistent with uncertain 
V opt (+4.7+15.0 kms -1 ). 

023.9— 02.3 (M 1-59) Single component with prominent broad wings. A CO spec- 
trum obtained with the JCMT 15 m telescope was previously reported by Gussie & Tay- 
lor (1995). Emission spike at V = — 8 kms -1 is IS. V Q consistent with uncertain V opt 
(+85.9 + 14.0 kms -1 ). 
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025.9—00.9 (Pe 1-14) Single component profile with series of narrow IS spike/dip fea- 
tures. V opt is unknown. 

025.9—10.9 (Na 2) Double-peaked profile. V Q consistent with uncertain V opt (+97.8 ± 
12.0 kms -1 ). 

031.7+01.7 (PC 20) Single-peaked profile. V opt is unknown. 

032.7+05.6 (K 3-4) Single peaked profile with broad wing on red side; IS spike/dip 
contamination on the blue wing from to —18 kms -1 . V opt is unknown. 

039.8+02.1 (K 3-17) Strong, remarkably broad, roughly flat-topped profile, possibly 
with additional weak wings; the dips in the profile at V = 13-18 kms -1 are IS. V opt is 
unknown. 

043.0— 03.0 (M4-14) Strong, asymmetric double-peaked profile, with broad wings. V Q 
consistent with uncertain V opt (+31.7+ 15.0 kms -1 ). 

047.1+04.1 (K 3-21) Tentative 5-a detection. V opt is unknown. 

055.3+02.7 (He 1-1) Tentative 5-cr detection. V Q differs by ~ 42 kms -1 from V opt ~ 
—35 kms -1 from Guerrero et al. (1999), but roughly coincides with one of the main optical 
components. 

055.6+02.1 (He 1-2) Flat-topped or concave profile with IS spike on the blue wing at 
V = — 8 kms -1 . V opt is unknown. 

056.0+02.0 (K 3-35) Flat-topped profile with IS spike on the blue wing at V = —8 kms -1 . 
A CO spectrum obtained with the NRAO 12 m telescope was previously reported by Dayal 
& Bieging (1996). V D agrees with H2O maser velocities in the core region (Miranda et al. 
2001). V opt (—9 ± 2 kms -1 , converted from LSR from Miranda et al. 2000) is offset by 
~ 12 kms -1 . 

060.1- 07.7 (NGC 6886) V Q consistent with very uncertain V opt (-35.8 ± 21.1 kms -1 ). 

060.5—00.3 (K 3-45) Strong, double-peaked profile; narrow dip at V = +15 kms -1 and 
spike at V = —32 kms -1 are IS. V opt is unknown. 

062.4—00.2 (M 2-48) Strong, (probably) double-peaked profile, with IS contamination 
of red side; narrow dip at V — — 11 kms -1 is also IS. V Q consistent with V opt ~ —15 kms -1 
in core from Lopez-Martin et al. (2002). 

067.9—00.2 (K 3-52) Double or single component with broad wings affected by strong 
IS contamination: dips in profile at V = —17, —4 and +4 kms -1 are IS. V opt is unknown. 
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068.3— 02.7 (He 2-459) Roughly parabolic profile with IS dips just outside the red wing. 
The very uncertain V opt (—72.0 ± 20.0 kms -1 ) differs from V Q by ~ 33 kms" 1 . CO(1-0) 
features at —9 and —15 kms -1 (+9 and +3 kms -1 LSR) reported by Zhang et al. (2000) 
are likely interstellar. 

069.6—03.9 (K 3-58) Strong, double-peaked profile. V opt is unknown. 

074.5+02. f (NGC 6881) Probably double peaked profile, but strong IS contamination 
(spikes) on blue side. V Q consistent with V opt (—14.4 ± 2.3 kms -1 ). 

079.6+05.8 (M 4-17) Flat-topped or concave profile. V Q consistent with very uncertain 
(-26.0 + 40.0 kms" 1 ). 

091.6-04.8 (K 3-84) Flat-topped profile. V opt is unknown. 

094.5— 00.8 (K 3-83) Asymmetric or double-peaked profile, affected by IS contamination 
on red side; dip at V — —57 kms -1 and just beyond the red wing are IS. V opt is unknown. 

104.4- 01.6 (M 2-53) Steep-sided, double-peaked profile. Deep IS dip at V = -43 kms -1 
affects the extreme red wing of the profile. V Q consistent with very uncertain V opt (—62.0 ± 
40.0 kms -1 ). 

107.6— 13.3 (Vy 2-3) Tentative 5a detection. Line is narrow for a PN but it could be 
the peak of a broader unseen component. It could be IS, but it is not seen in any of the 
surrounding 8 pixels, and V Q consistent with well determined V opt (—49.5 ± 3.8 kms -1 ). 

119.6— 06.7 (Hu 1-1) Single-peaked, asymmetric profile. The well determined V opt 
(—53.7 ± 3.0) is offset by ~ 14 kms -1 from the V Q peak, in the direction of the extended 
wing. 

130.4+03.1 (K 3-92) Double-peaked profile. V Q consistent with V opt (-61.7+2.8 kms -1 ). 

149.0+04.4 (K 4-47) Strong, roughly parabolic profile. V opt consistent with V opt ~ 
—37 kms -1 from Corradi et al. (2000). 

153.7- 01.4 (K 3-65) Double-peaked profile; spike at V = -26 kms -1 is IS. V opt is 
unknown. 

184.6+00.6 (K 3-70) Flat-topped profile. V consistent with V opt (+26.9 + 4.5 kms -1 ). 

359.8+06.9 (M 3-37) Strong, double-peaked profile; spike at V = — 5 kms -1 is IS. V D 
consistent with V opt (—74.2 ± 10.2 kms -1 ). 
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5. 



Discussion 



5.1. 



Detection Rates 



The observations reported here significantly increase the available data on the CO en- 
velopes of PNe. The compilation by Huggins et al. (1996) lists 44 PNe detected in CO, and 
only a handful of other PNe have been detected in searches in the intervening period (e.g., 
Josselin & Bachiller 2000). One focus of the earlier work was large, highly evolved PNe with 
low surface brightness. The new sample has a higher proportion (85%) of PNe with moder- 
ate to high surface brightness (Tb > 65 K using the benchmark of Section 2.1), compared 
to 64% of the earlier compilation, and roughly doubles the number of relatively young PNe 
with known CO envelopes. The 70 sensitive upper limits reported here also improve and 
extend the earlier data on PNe where the molecular component has been destroyed. 

The overall detection rate of the new sample is 36%, and the rate shows no significant 
variation with Tb. On the other hand, the detection rates for the spectroscopic groups 1-4 are 
highly non-uniform at 15%, 71%, 36%, and 18%, respectively; this variation is significantly 
different from that expected on the hypothesis of a uniform distribution at the 0.1% level. 
These results are consistent with the view that large intrinsic variations of the CO properties 
of the PNe (Huggins et al. 1996) dominate over evolutionary effects for the ensemble as a 
whole. They also indicate that the presence of a CO envelope is correlated with higher [N II] 
emission in the nebula. This is likely to be mainly a direct effect in which [N II] is enhanced 
in dense interface gas, or shocked gas often associated with bipolar flows and dense tori. It 
may also be partly indirect, in that higher mass progenitors, which often have an enhanced 
N abundance, preferentially have more massive molecular envelopes. These will be discussed 
in more detail in a separate paper. 



The CO spectra of the 40 detected PNe provide a useful and homogeneous set of PN 
radial velocities (V Q in Table 1). In the cases where the CO spectrum is relatively clear of IS 
line contamination and the S/N is good, the uncertainty in the measured velocities is small 
~ 1-2 kms -1 , and is limited by the intrinsic asymmetry in the line profiles. Even in cases 
with severe IS contamination or low S/N, the uncertainty is almost always < 10 kms -1 . 

High or medium quality (uncertainties < 10 kms -1 ) optical measurements of the radial 
velocities of the detected PNe are available in only 18 cases (see Section 4). For this set, 
the mean value of V Q — V ovt = 5 kms -1 with an rms scatter of 12 kms -1 . Thus the CO 



5.2. 



Radial Velocities 
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and optical measurements provide similar, robust estimates of the systemic velocity, with 
part of the residual scatter being real, reflecting differences in the structure of the ionized 
and molecular gas in the PNe. In the other 22 cases where CO is detected but the optical 
uncertainties are large, > 10 kms -1 (10 PNe), or are non-existent (12 PNe), the CO data 
provide the preferred or the only available radial velocity measurements of the PNe. 

5.3. Line Profiles and Widths 

A striking feature of the CO spectra shown in Fig. 1 is the variety of the line profiles. 
They can be broadly classified into 4 basic types: (I) parabolic, similar to the CO spectra 
typically seen in the precursor AGB envelopes; (II) flat-topped; (III) double-peaked; and (IV) 
triangular, usually asymmetrical. There are also intermediate types: e.g., with a distinct but 
not deep concavity, between types II and III; and with strongly asymmetrical double-peaks, 
between types III and IV. 

These different profiles probably arise from a combination of effects. The first is the 
evolution of the structure of the envelope from the relatively complete, spherically symmetric 
distribution seen around AGB stars, to one in which the envelope has been hollowed out by 
dissociating and ionizing radiation and/or shaped by the action of jet outflows (see Section 
1) into bi-cones, rings, tori, etc., which may be oriented at different angles to the line of 
sight. Red-blue asymmetries in the profiles then arise from back- front asymmetries in the 
distribution of molecular gas. A second effect is the likely decrease in the opacity of the CO 
line as the CO is destroyed. A third effect is the size of the PN cavity relative to the telescope 
beam size. If the third effect were dominant, one would expect that the profile shapes would 
be correlated with the angular diameter of the PNe. 068.3—02.1 is among the smallest PNe 
(1.3" Aaquist & Kwok 1990), and is the case with the most parabolic profile, consistent 
with the emission arising in a relatively undisturbed envelope. However, there is no general 
correlation between the profile shapes and the angular sizes of the ionized nebulae for the 
sample as a whole. The prevalence of asymmetries, and strongly double-peaked profiles 
even among small nebulae, suggest that the structural effects in the envelopes are of major 
importance. 

A second, striking feature of the CO spectra is the large widths of the lines. The 
distribution of AV listed in Table 1 is shown in Fig. 2, where the lower limits for AV 
obtained from Gaussian fits are included. For comparison, Fig. 2 also shows the distribution 
of CO line widths (twice the expansion velocity) of the sample of 65 bright C-rich AGB stars 
reported by Olofsson et al. (1993); the distribution of the widths of OH/IR stars is similar 
as discussed in that paper. The difference between PNe and AGB line widths is clear from 
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the figure. High velocity gas is rarely seen in AGB envelopes (median width = 25 kms 1 ), 
but is the norm in PNe envelopes (median width > 45 kms -1 ). 

The expansion velocity of the ionized gas in PNe is also larger than the expansion 
velocity of AGB envelopes, typically by a factor of about two. The mean expansion velocities 
given for PNe in the catalog of Weinberger (1989) are 20.3 kms -1 for [O III] and 22.5 kms -1 
for [N II], and the corresponding line widths fall near the center of the distribution of CO 
line widths shown in Fig. 2. In PNe the large expansion velocities of the ionized gas are 
usually ascribed to acceleration by the pressure in the nebulae caused by photo-ionization 
or fast winds. In the case of the CO lines, the large line widths of the PN spectra in Fig. 1 
are mostly formed by extended line wings, e.g., M 1-59 (023.9—02.3). Similar wings are 
commonly seen in proto-PNe before the onset of ionization, and in the few cases of young 
PNe where the CO line wings have been studied at high angular resolution, e.g., M 1-16 
(Huggins et al. 2000) and He 3-1475 (Huggins et al. 2004), it is found that the wings arise 
from molecular gas entrained in the sides of high velocity jets. Thus although the kinematics 
of the molecular gas may be affected by thermal pressure in the ionized gas, the frequency 
with which wings are seen in the CO spectra reported here probably reflects the frequency 
of jets in PN formation. 

In a few of the cases reported here, the whole CO line appears to be broadened. The 
most extreme case is K 3-17 (039.8+02.1) which shows a steep-sided, roughly flat-topped 
profile (with some IS contamination near the center) that is 90 kms -1 wide. This profile is 
unique among PNe and the object deserves further study. 

6. Conclusions 

This paper reports a sensitive survey of CO emission in PNe. 110 PNe were observed and 
40 were detected. The success of the survey demonstrates the usefulness of array receivers 
for observing CO in PNe or other small sources lying close to the Galactic plane. The results 
substantially extend the available data on the molecular gas in PNe for statistical study, and 
provide individual cases for further detailed observation. 

We thank the staff of the IRAM 30 m telescope for help with the observations. This 
work has been supported in part by NSF grants AST 99-86159 and AST 03-07277 (PJH), 
and by the Spanish MEC under grants AYA2003-7584 (RB), ESP2003-04957 and FEDER 
grants (PP). 
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Fig. 1. — CO J = 2 — 1 spectra of the PNe detected in the survey. The PNe are identified 
by the galactic name in each panel and are ordered by Galactic longitude. See section 4 for 
comments on individual spectra. 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 2. — Line width comparison of PNe and AGB stars. Lower panel: Distribution of the 
CO line widths of the PNe detected in the survey. The darker regions denote lower limits to 
the CO widths. Upper panel: Distribution of the CO line widths of 65 bright carbon stars. 
Data from Olofsson et al. (1993). 
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Table 1. CO Survey of Young Planetary Nebulae 
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Table 1 — Continued 
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Note. - Column [9]: d, detection; n, non-detection; t, tentative detection; gl-4, spectroscopic 
group; s, severe interstellar contamination 



